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We obtain modified dispersion relations by requiring the vanishing of determinant of inverse 
of modified photon propagators in Lorentz invariance violation (LIV) theory. Inspired by these 
dispersion relations, we give a more general dispersion relation with less assumption and apply it to 
the recent observed gamma-ray burst GRB090510 to extract various constraints on LIV parameters. 
We find that the constraint on quantum gravity mass is slightly larger than the Planck mass but is 
consistent with other recent observations, so the corresponding LIV coefficient £1 has reached the 
natural order (0(1)) as one expects. From our analysis, the linear LIV corrections to photon group 
velocity might be not excluded yet. 
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I. INTRODUCTION 

c : 

a : , 

)—} , Lorentz invariance violation (LIV or LV) has been intensively investigated both theoretically and experimentally 
in recent years. The revival passion of relativity violation in theoretical construction originates from the attempt to 
t-H | compromise general relativity with quantum mechanics. On the other hand, the experimental searches may provide 
us with concrete evidence to sift a most hopeful candidate of quantum gravity from a vast number of theories. 

From theoretical aspect, some theories expect LIV to happen at high energies. For example, spontaneous Lorentz 
symmetry breaking may happen in string theory as the perturbative string vacua is unstable, thus some tensor fields 
generate nonzero vacuum expectation values [l|. The breaking of Lorentz symmetry also happens in other frameworks, 
such as loop gravity 0], foamy structure of spacetime 0], torsion in general gravity etc.. More recently, Hofava 
proposed a power counting renormalizable theory of gravity |5| with a "dynamical critical exponent" z to characterize 
the anisotropic scaling properties between space and time. While Lorentz symmetry is breaking at high energies, it 
restores when this dynamical critical exponent flows to z = 1 at low energies. There are also some other proposals, 
such as the so called double special relativity Q, which preserves relativity principle with a nonlinear realization 
of Lorentz group, thus conventional Lorentz symmetry is also broken. One striking consequence of LIV is that the 
' photon propagation speed is no longer a unique constant, generally, it depends on energy and propagation direction. 
, These theoretical investigations have promoted various experiments to search for the deviation from conventional 
■ linear dispersion relation for photons [7]. However, as the possible violation effects for photons must be very tiny, the 
\ detection of these effects present a significant challenge to experimentalists. In addition to improve the precision of 
measurements to find any possible evidence of LIV, we should also take efforts on searching for certain accumulating 
^\ ' processes to amplify these tiny effects. Such idea has already been proposed on the observation of certain astronomical 
objects such as gamma-ray bursts (GRB) [!, H|, pulsars (9[ and active galactic nuclei (AGN) [Io| . etc., and the tiny 
LIV effect could manifest itself through the observation of rotation of linear polarization (birefringence) [ll[ or time 
• t-H . of flight lag [1] for photons with different energies. 

The paper is organized as follows. In Section 2, we review certain modified photon dispersion relations derived from 
several LIV models, including standard model extension (SME) with power counting renormalizable operators [TTj j . 
effective field theory with dimension 5 operators [l2| and Hofava's anisotropic U(l) theory. In Section 3, we focus on 
time of flight analysis of GRB and try to extract some LIV parameters from the recent observation of GRB090510 13]. 
We briefly discuss the time of flight analysis of photons from cosmological distant objects, then we give a general 
dispersion relation used conventionally in the astrophysical analysis of LIV [T^. This general dispersion relation 
contains those terms derived from the models in Section 2 as special cases. We then extract constraints to linear LIV 
parameters from GRB090510 to O(0. 1), improved by 1 or 2 order of magnitude than those in [3] and [T(J[l5j]. From 
the analysis of the time-lag formula we find that it is hard to significantly improve the constraints from this simple 
and rough analysis, unless other time-lag effects (like source effect [lq . which is a major uncertainty in the time of 
flight analysis) can be clarified or other methods will be used. 
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II. PHOTON DISPERSION RELATIONS 
A. Background tensor field induced LIV 

A systematical treatment of LIV to incorporate particle standard model with power counting renormalizable La- 
grangian, called standard model extension (SME), was proposed by Kostelecky and Colladay in Ref. [llj, where the 
photon sector reads 

£ P hoto„ = ~\V V - \(k F ),x^F KX F^ + ^(k AF ) K e KX ^A x F^. (1) 
From ([T]), we deduce the equation of motion below 

d a F m + (k F )^ aPl d a F^ + {k AF ) a e^ aM F^ = 0. (2) 
By expressing @ in terms of 4-vector potential and assuming that the Fourier decomposition 

A^(x) = a f _ l (p)exp(~i p ■ x) (3) 
is still reliable, we express ([2]) in the momentum space as 

AV(p)a» = 0, (4) 

where 

M^(p) = V^v P 2 -P^Pv- 2(k F ) flKX vP K P X - 2i(k AF ) K 'e i _ lK xuP X - (5) 

We impose gauge fixing condition and require the determinant of the reduced matrix to vanish, then we can obtain 
an implicit function p°(p), which is an eighth order-polynomial in p a . Otherwise one can verify that the determinant 
of M^y vanishes as a consequence of gauge invariance of equation ([5]). For example, we use Lorentz gauge 

d a A a = 0, (6) 

in momentum space, i.e., 

Pa a a (p)=0. (7) 

So we have the gauge fixed reduced matrix 

Mf u (p) = p 2 - 2(k F )^ x „p K p x - 2i{k AF ) K e^p x . (8) 

For our purpose, we just try to extract a simplified result by assuming that (particle) rotational invariance still holds 
regardless of the explicit violation of Lorentz symmetry, i.e., only (k AF )° and a (a combination of {k F ) K \^ v , for 
details, see Appendix or [13]) are nonzero. With this assumption, we have the following matrix: 

M rcd {p) = 



9-9 ft 1 9 n ^ 

p — ap ap p ap p ap p 

ap°p 1 -{p 2 + a((p ) 2 +f -{p 1 ) 2 )) ap 1 p 2 + 2ik° AF p 3 ap V - 2ik° AF p 2 

ap°p 2 ap 1 p 2 - 2ik AF p 3 -{p 2 + a((p ) 2 + p 2 - (p 2 ) 2 )) ap 2 p 3 + 2ik AF p 1 

ap°p 3 ap 1 p 3 + 2ik AF p 2 ap 2 p 3 - 2ik AF p 1 -(p 2 + a(p ) 2 + p 2 - {p 3 ) 2 )) 



(9) 



Its determinant reads 

det(M rcd (p)) = {i(k AF ff - ((1 + «)(P°) 2 - (1 - a)f) 2 } (1 + a){p 2 ) 2 . (10) 

By requiring det(M red (p)) = (otherwise there would be no solution for a photon field), we have two dispersion 
relations, one is the conventional p 2 = and the other is 

(P°) 2 = Tf^y ((! - <*)? ± Maf\pA) ■ (11) 
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We can also use another equivalent method to obtain these two dispersion relations. First, we rewrite (pj in an 
explicitly quadratic form in the photon field, i.e., 

£ P hoto„ = - \f^ v - - A {k F ) KX ^F^F^ + \{k AF )^ Xp »A x F pv - i(0 • A) 2 
= -\ d^Av [F^ + 2(k F ) KX ^d K A x + ■ A - 2e KX ^(k AF ) K A x \ 



= total derivative + - A u (D F L y x A x , 
where we have added gauge fixing terms — j^{d ■ A) 2 in (fT2"j) and 



(D 



~ 1 ) vX = (ttf x - d u d x (l - |) - 2(k F ) upnX d p d K - 2e^ KX (k AF )^d K ^j 



Using the same Ansatz (|3]), we define a matrix £ in momentum space 



P~T)v P - (1 - ^)PuP P ) + 2{k F ) vlllip p t "p K - 2ie vpKp p tl {k AF y 



(12) 



(13) 



(14) 



From the conventional free field theory, the differential operator inside the two fields in the quadratic form of certain 
Lagrangian (e.g., (fl~2"|) ') is just the inverse of free field propagator in position space (see [TH or [l^| ), thus (fl4"|) is just 
the inverse of photon propagator expressed in momentum space. We know that generally the inverse of propagator 
is just the dispersion relation, from which one can find the pole of the corresponding particle, so we expect that 
the determinant of (I14[) in case of £ — > oo (i.e., without gauge fixing) is zero. Similarly, we can find the explicit 
dispersion relation in a special gauge by choosing the corresponding specific value of £. For example, we find that 
£(p)„ p |£ = i = — M ei {p) lJp (£ = 1 is just the Lorentz gauge (JSJ) used to obtain ©, and this choice can avoid the 
inequivalent gauge choice comparison. As pointed out in [llfl , different gauge choices are inequivalent with each other 
in the LIV electrodynamics). So in the rotational invariant case, this matrix can also lead to (|11[) and the conventional 
dispersion relation. We mention here that similar method to obtain photon propagator in the SME framework has 
also been obtained recently in [2(J, with a more systematic and complete treatment. 

The leading order nonrenormalizable LIV operators (dimension 5) were systematically studied in [l2| , where Myers 
and Pospelov also introduced explicitly a timelike four- vector n a to take LIV into account, thus this theory can be 
regarded as a leading nonrenomalizable part of SME. Since we are only interested in the study of the consequence of 
LIV to the propagation of GRB, we focus our attention only on photon field there. The corresponding Lagrangian is 



SC 



photon 



2M P , 



e^ Kp n a F ap n-d(n K F p „). 



We write it in another equivalent form, i.e. 



-■photon 



- A v Urf" - 



2£ , 
Mpx 



d(n ■ d n K d„e^ Kp + n K d p d a t vpKa n p ) A p + total derivative 



(15) 



(16) 



where we have added the Lorentz gauge fixing term. Then by performing the same procedure as before, we have the 
reduced inverse of propagator 



2i£ 
~Mp~i 



(e^p 2 oPfl + e^ 0a p oPp p a S p ) 



(17) 



when expressing explicitly the time-like four-vector n in a preferred frames as n p = (1,0,0,0). Then by imposing 







detn(p) = det 



f-p 2 

p z 

o ijg-OW p 2 

we obtain the dispersion relation 

(p°r=f± 



\ 



JL 

'M F 



{p0) 2 p l 



P 4 [(^P) 2 (P°) 4 -P' 



M 



pi 



0, 



(18) 



2£ 



(P°) 2 \P\, 



(19) 



which was obtained in |12j | plus the conventional one p 2 = 0. 
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B. Anisotropic scaling induced LIV 

Now wc turn to another framework of LIV proposed recently by Horava Q . His original proposal was to provide a 
UV completion of quantum theory of gravity. Lorentz symmetry appears naturally in this theory when the dynamical 
critical exponent flows to z — 1 at low energies. While at high energies, space and time present anisotropic scaling 

t^XH, r^Xr, (20) 

thus Lorentz symmetry breaks down. However, this formalism does not break spatial isotropy, thus there is no need to 
assume a special background field configuration to realize rotational invariance, unlike the background tensor formalism 
discussed above. Aside from gravity, Hofava also constructed an anisotropic Yang-Mills theory with critical spatial 
dimension D = 4 [f| . As Chen and Huang recently gave a general construction of bosonic field theory demonstrating 
this anisotropic scaling behavior [21] . we follow this new approach instead of [5|. In the new formalism, the photon 
action reads 

S=lf dtd D x\ ( E 2 ~Y J-fVl)" w^'" n ; d 2n *F J ). (21) 

For simplicity, we consider the case z — 2 and D — 3. Then one immediately reads from the action that the scaling 
dimensions of the couplings are 

[g E ] s = \{z-D) + l, [Xj, n ] s =z + D+~(z-D-2)J-2n. (22) 

Thus in this case the renormalizable condition ([<?,e] s > : z > D — 2) for E is automatically satisfied. Actually, it is 
superrenormalizable. If the critical dimension (i.e., [gE\s = 0) is D — 3, then z must equal to 1, which just corresponds 
to the conventional Lorentz invariant gauge theory. Renormalizability also imposes the condition [Aj „] s > 0, and for 
a free field theory, J = 2, n < z — 1 = 1. For simplicity, we set A2.0 = 5, then the free Lagrangian (with gauge fixing 
term) is 

£frcc - 4r (e 2 - ^FijF** - ^(diFik ■ djF jk + diFjk ■ dtFjk^j - ^(d- A) 2 



±^A V I (of" - d v d»{l - i)) - ^A(A5 fcj - dkdiW^ A p . (23) 



By performing the same trick, we can obtain 



T( P r = -p 2 ^" + -^8l5 P 3 { PkPj - fS jk ), (24) 



and the corresponding dispersion relations read 



P 2 =0, { P y=f{l + ^f). (25) 



III. TIME OF FLIGHT ANALYSIS OF GRB IN LIV THEORY 



In this section, we discuss the LIV effect on the observed GRBs, focusing especially on the time of flight of 7 rays. 
Gamma-ray bursts (GRBs) are sudden, intense flashes originating from distant galaxies with cosmological distances, 
and they are the most luminous electromagnetic events we ever known. As already shown in the above formulas, 
LIV can modify conventional Maxwell equations and hence leads to modified dispersion relation in addition to the 
conventional one. However, due to the large mass scale suppression, these LIV effects must be very tiny to account 
for the conventional stringent terrestrial test. Fortunately, as first pointed out in [3j, the cosmological origin plus high 
energy and the millisecond time structure of GRB make GRB an ideal object to observe the possible minuscule effects of 
LIV. Actually, many known stringent constraints to LIV parameters were drawn from astronomical observations, such 
as AGN ultrahigh-energy cosmic rays (UHECR) 22], CMB [H,[2l|, etc.. The LIV induced modified dispersion 
relation can lead to many interesting phenomena. The most apparent consequence is the frequency dependence of 
photon group velocity, though this is not always the case. For example, if only fci? 7^ in the SME framework, photons 
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propagate independently with their energies. So if photons with different energies are emitted simultaneously, this 
frequency dispersion of group velocity then leads to the so called time-lag phenomenon. In addition to time lag, 
certain models, e.g. SME, indicate that photons with independent polarizations obey distinct dispersion relations. 
This was demonstrated in the above two models, see (fTTj) and ([T9")) . All these models involve a conventional mode 
with an extraordinary helicity dependent one, thus can lead to the so called vacuum birefringence effects [111 |26| . The 
tiny changes in polarization grow linearly with propagation distance and hence can be accumulated to be observable 
for cosmological sources. This can provide a sensitive probe to LIV [Til [23Tj26| . Aside from purely kinematic effects, 
the tiny LIV correction to dispersion relation can also dramatically change the thresholds of high-energy particle 
reactions, hence leads to distinct particle spectrum of UHECR with respect to that of Lorentz invariance cases. The 
observation of this spectrum can provide a unique signature of LIV [27j | . On the contrary, the nonobservation of these 
effects can put very stringent constraints to LIV parameters [22j|. Below we will primarily discuss the GRB time-lag 
caused by LIV. 

First we make a brief review of the formula used in the description of GRB photon time-lag with respect to the 
source redshift. For an isotropic and homogeneous universe, one can derive a differential relation 

dt = dz (26) 

h {\ + z)^jn h + n K (i + z) 2 + n M (i + z y + n R {\ + z y 

from the Friedman equation 



a ?+ K 
a a 6 



see [28[ for details, where the present-day Hubble constant Hq ~ 71 km/s/Mpc and Qk = for a nearly flat universe. 
The matter density Qm — 0.27, radiation density Qr ~ and vacuum energy density fl\ ~ 0.73 are the cosmological 
parameters evaluated today. We note that ([26]) is the standard result derived from general relativity which is a locally 
Lorentz invariant (LI) theory. So in dealing with photon time-lag below, we implicitly assume that the gravity side 
is untouched. Though a unified treatment should also include the change of gravity due to possible LIV effects hence 
may also change ([26]) and the time-lag formula used below. 

Then we make a general assumption of photon dispersion relation 

E 2 =/(p;M,{&}), (28) 

where f(p;M, {^}) is a general function of p (p — |p|), some unknown large scale M relevant to LIV and a set of 
parameters {^}. Inspired by the dispersion relations ([19]) and ([25]) derived from particular models discussed above 
and the fact that LIV corrections must be very tiny at low energies, we assume that the expansion of ([28]) around the 
conventional dispersion relation E 2 = p 2 is 

N 

£ 2 = P 2 (i+]>>(^r), (29) 

1=1 

where N is a large number marking the precision of our expansion. Thus (| L9[) and (1251) can be regarded as just two 
special cases of ([29]) . i.e., only ^ ^ and ^ ^0 respectively (where ([T9]) just adds helicity dependence assumption 
of LIV) . In addition to its generality, the reason for beginning with ([2"9"]l instead of those particular models is that 
various experiments have already ruled out (TT9l to a convincing level (see [29[ and a recent review [30[). Indeed, 
dimension 5 LIV operators have suffered very stringent constraints both from frequency-dependent birefringence test 
with GRB pnl |32T]. Crab Nebula [33[, and the UHECR spectrum analysis under special assumptions [22} (i.e. LIV 
corrections to electron dispersion relations are smaller than those of photon ones, e.g., in the Liouville string models 
of foamy structure of space-time [Hj] , where only neutral gauge bosons receive quantum-gravity corrections) . 

By taking into account of the expansion of universe [35| and the assumption that gravity side remains intact, the 
time-lag led by modified dispersion relation ([2"9")) with leading order LIV correction of order n is 

St=^J-^ r^^dz', (30) 
2 ^ n M n J a h(z') ' v ' 

where 



h(z) = H ^Q A + Q K (l + zf + fiAf (1 + zf + n R (l + z)\ 



(31) 
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Eq is the redshifted photon energy observed on earth. SEq = EJ l — E h l , where Ei and Eh denote lower and higher 
energies of observed photons respectively in the time delay. The n-th order correction corresponds to dimension n + 4 
LIV operators. This can be seen from another formula 

St = Swd ~* f {1+ h(zV dz ' E oY jm (n)k[% m , (32) 



which is suitable to the analysis of time-lag in the SME framework given recently by Kostelecky and Mewes [36j j . 

Before we discuss the linear (n = 1) and the quadratic {n — 2) corrections to the photon dispersion relation, we 
first utilize the observed time delay in GRB090510 located at redshift z — 0.903 ± 0.003 to give a rough estimate to 
photon mass. As is well known, photon mass is represented by dimension 2 op erator (quadratic in photon fields) and 
may spoil gauge invariance (this is not the case in Chern-Simons theory [19| in space-time dimension 3). However, 
the presence of photon mass does not necessarily implies LIV. In a LI theory (Proca's theory), the existence of a 
unique speed c could be regarded as the limiting speed of light for arbitrary-high energy photons. So photon mass 
is absent in a LIV theory (e.g. SME J3jj ) if gauge invariance is still valid. But we can still give a rough estimate of 
its magnitude by using (|32|) and the fact that the bulk of the photons above 30 MeV arrived 258 ± 34 ms later than 
those below 1 MeV in the observation of GRB090510 : 

oY jm (n)kff )jm < 1.4801 * 10~ 24 GeV 2 . (33) 

jm 

This could be translated to the photon mass bound as m 7 < 1.217 * 10~ 3 eV, much larger than the mass upper 



bound given in |37|, m 7 < 1 * 10 -18 eV. This confirms the remarks given in [38| : "(departures of electrostatic 
and magnetostatic fields from the gauge invariant one) give more sensitive ways to detect a photon mass than 
the observation of velocity dispersion." Of course, one can obtain an effective mass bound comparable to this as 
m 1 < 1.217 * 10 -19 eV [24|, but the origin is different. From the bound derived we see that one would need a much 
larger photon mass to explain the time of flight data if one does not introduce the LIV effect (or other effect, e.g. source 
effect). More than that, as the presence of photon mass indicates that high energy photons propagate faster than low 
energy ones, this time advance of high energy photons may cancel possible time-lag induced by certain LIV models 
(£i < in (|29j) . see also ([35])). Thus mass effects may conspire with LIV effects to produce a nearly nonobservation 
of time-lag in certain time of flight analysis [39j . On the contrary, in some scenario with £j > 0, the time-lag might 
be caused by the combined effects of mass and LIV, thus the situation is still complicated. Fortunately, due to the 
high precision laboratory experiment [40| (constrain m 7 to 10 -17 eV level), those scenarios mentioned above do not 
happen and we can safely ignore the mass effects in our discussion about LIV constraints drawn from the time of 
flight analysis of GRB090510. 

Without the trouble of possible mass effects, we can then securely discuss LIV effects in the time-lag phenomena 
below. As a byproduct of (|3"2"|) . we give a rough estimate to mass dimension 3 LIV operators 



oY jm (fi)k$ Jm < 1.1558 * 10- 21 GeV. (34) 

jm 

We see that this bound is comparable to that obtained from the LIV effects on Schumann resonances in a natural 
earth- ionosphere cavity [41[ , though it is much weaker than the other astronomical constraints |23l [24J (which are 
constrained to less than 10 -43 GeV). 

Then we turn to nonrenormalizablc LIV operators but using formula (|30|) of time-lag instead, as it is more suitable to 
our simple analysis based on general dispersion relation (|29p . As usual, we only discuss energy dependent corrections 
to photon group velocity to the quadratic level. Before looking into details, from (|29l) we derive modified group 
velocity 

9 dp E \ 2 y ' M 
_ + + 2) &(£)«) 

^l + ^ + l)^) 4 , (35) 

where Einstein sum over index i from 1 to N is indicated. Then by the same procedure in [35| . we can give a time-lag 
formula which is accurate to 2nd order of and first order in Az as 

Ei-Ei f (1+^0 3 ^ Ef-El r z (l + z')\ , 

61 = J ~WY + s (46 - Cl) ^ 2_ 1 ~WV dz ■ (36) 
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It can be seen that (|36[) is consistent with (|30l) when £2 = and £1 = respectively for linear (though the linear 
correction is calculated to second order in the large mass suppression, it will be checked that this can not improve the 
linear constraints any more, thus in practical calculation, (|30|) is enough) and quadratic corrections to group velocity. 
For linear energy dependent correction to 1st order, we derive from the most conservative claim that, the bulk of the 

photons above 30 MeV arrived 258±34 ms later than those below 1 MeV pj], the LIV scale ^ 5.02689*10 16 GeV, 

which is 3-order less than the Planck scale if |£i| is of order 1. However, if utilizing the more stringent claim that, a 
single highest detected photon from GRB090510 with 31 GeV arrives 0.179 s later than the main LAT emission above 
100 MeV, we can deduce a significant higher quantum gravity mass scale 



M 



-7.72017 * 10 19 GeV, 



(37) 



where the minus sign indicates the fact that photons with higher energies propagate slower than lower ones as already 
mentioned in the discussion of photon mass (£1 < 0). By direct calculation of solving 2nd order equation of |^ (i.e. 

setting £2 = in (J36J) ) , we find that this can not improve the result any more as mentioned. This illustrates that in a 



rough estimate of linear correction to group velocity, there is no need to take into account 2nd order correction of ( |^) 2 
as (|36|) . The result ([37]) means that linear correction gives a LIV mass scale ne arly 6.32 Mpi if £1 ~ 0(1), which is very 
close to that of [llj]. Of course, if one chooses other data from the Table 2 in [13], one can obtain the same conclusion 
that quantum-gravity mass scale is significantly above the Planck mass (at most of order 102 Mpi [l3j]) from this simple 
analysis. This conclusion is nothing more than a translation of the claim that the constraint on the linear energy 
dependent LIV parameter £1 can be placed in the range 10 _1 ~ 10~ 2 , if we regard M ~ Mp\. If this constraint 
is confirmed by other astrophysical observations, then it puts the constraints at least 2 order of magnitude stronger 
than [l(| and [15j], which gives |£i| < 17 and |£i| < 58 respectively. However, these constraints are not stronger enough 
as those obtained in [22j extracted from the UHECR spectrum and those in [3l[ from the frequency-dependent helicity 
observations of GRB930131 and GRB960924. However, we note that those most stringent constraints (£1 < 10~ 14 ) 
up to now rely either on particular assumptions (see [22j ) or helicity dependent models, e.g. (|19[) . Thus it is necessary 
to put the constraints obtained from (f2l))) on the linear LIV parameter £1 to the same level (still a hard task as a span 
of 12 orders to be conquered) from future observations. If so, we can finally make a more firm claim that dimension 5 
LIV operators can be excluded firmly (30} . Then we may reach the conclusion in the near future that either Lorentz 
symmetry is exact, or at high energies there are some other symmetries such as SUSY plus CPT to protect our low 
energy theory from receiving CPT odd corrections [29| . 

For quadratic energy dependent correction, i.e. £1 = 0, we obtain from the most conservative claim mentioned 
above the constraints on quantum-gravity mass scale M ~ 5.84718 * 10 7 GeV if £2 ~ 0(1). While for the single 

31 GeV event, we obtain the constraint as ^ 5.26767 * 10 21 GeV 2 , thus M ~ 7.25787* 10 10 GeV. It is obvious 

that the constraints obtained from quadratic correction are much weaker than those from linear one as quadratic 
correction being suppressed more than one power of M. Thus in the future we should take more efforts to the 
search of more stringent constraints on quadratic LIV correction to photon group velocity. As mentioned above, 
the quadratic correction is produced by dimension 6 operators, which are the leading order nonrenormalizable CPT 
even LIV operators. As we known, various current constraints to dimension 6 operators are also much weaker than 
dimension 5 ones |3l| . 

Before we close this section, we observe that our results are similar to that obtained recently in pl| and pj|j EH . 
We give these results in the table below: 



Source 


Mkn501 [15] 


PKS 2155 - 304 [10] 


GRB080916C [14] 


GRB090510 [13] 


redshift 


0.034 


0.116 


4.35 


0.900 


St(a) 


240 


27 


16.54 


0.179 


E h to Ei{GeV) 


10 4 to 250 


600 to 210 


13.22 to 10" 3 


31 to 0.1 


M (GeV) 


6.06 * 10 17 


7.51 * 10 17 


1.55 * 10 18 


7.72 * 10 19 


t total 

at 


6.01 * 10 13 


1.72 * 10 15 


2.34 * 10 16 


1.29 * 10 18 



TABLE I: Where the first three rows below the source row are the data given from [TcL Il3l - fl5| respectively, the last 2 rows are 
the linear quantum-gravity masses and total time to time-lag ratios calculated from (|30|) . 



We find that the rough estimates about the linear mass scale are consistent with those given by the reference above 
and the order of magnitude of the linear mass scale ranges from 10 17 to 10 19 . It is easily seen from Table I that 
one can approach this large magnitude with the advantages both from the quotient of total time to time-lag (which 



8 



originates from the cosmological distance and short pulse nature of GRB) and the large absolute energy difference 
(range from GeV to TeV). To further constrain the linear order quantum- gravity mass hence the coefficient ^ in 
the future, we need to amplify the large ratio of total time to time lag, as photons with energies much higher than 
already observed (TeV) can not reach us from cosmological distance due to the pair creation interaction with infrared 
background photons. One way to amplify the large time ratio is to exclude other non-LIV induced time-lag factors, like 
different response time of detectors [42j (which slightly increase linear quantum-gravity mass scale obtained in [13]), 
or one turns attention to other methods like spectrum analysis (22j . otherwise the constraints can not be improved 
significantly. Further more, a statistic analysis by taking into account of statistic error [15J and a multi-source analysis 
to calculate the correlation between distance and time-lag [8j will make the results more concrete. 



IV. CONCLUSION 



In this paper, we reviewed several modified dispersion relations from standard model extension and Hofava theory in 
the photon sector. Dispersion relations are derived consistently from the inverse of photon free propagators, without 
taking quantum corrections into account. Inspired by these dispersion relations we give a more general one (|29[) 
to avoid some particular assumptions (e.g. helicity dependence). Then we apply this relation to the time of flight 
analysis of recently reported GRB0905I0. We obtain constraints on the linear LIV energy dependent coefficients to 
the level of £1 ~ O(.l), which is equivalent to the statement that the relevant linear quantum-gravity mass scale 
is M ~ 7.72 * 10 19 GeV. Using the same method we also get the quadratic mass scale M q ~ 7.26 * 10 10 GeV, 
denoting a much loose constraints to dimension 6 operators. These results are consistent with those in fiot Il3l - [l5| . 
As a byproduct of the time-lag formula (|32[) , we point out that one can safely ignore the photon mass effect in the 
discussion of LIV effects in the time-lag analysis of GRBs due to the stringent terrestrial constraints on photon mass, 
and we obtain a constraints J2j m oYj m (n)k^- m < 1.1558 * 1CT 21 GeV to the dimension 3 operator. 

From our analysis, we find that though the constraints obtained are far from reaching those from the spectrum 
analysis [22j and those from the helicity dependent analysis [3lj , our analysis relies little on extra assumptions except 
the expansion (|29j) and the formula (|30[) . Thus to exclude the linear order quantum- gravity correction to photon group 
velocity is still too early as long as the constraints to the general linear order correction (|2"9")l have not approached the 
same level as in [3(| ■ We find that the results have already reached the precision of probing Planck mass scale or even 
higher, slightly better than [13, [l5j]. If one can largely clarify other time lag uncertainties like (42j, the constraints 
could be improved more. 
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Appendix 

The previous referred parameter a, is just one of the parameters defined from various combinations of {kp) K \^ u . 
These definitions arise for convenience from the consideration of the symmetry of this tensor. From the Lagrangian 

SL = ~(k F ) K x„ v F KX F^, (38) 



we find that (k F ) K \^ lJ is antisymmetric to the two indices k\ and \j,v respectively, and is symmetric to the interchange 
of these two pairs of indices. As we do not want to include a conceivable 0-type term proportional to ^e K x t _ lu F KX F^ l/ , 
we require that e ti x^(k F ) KX ^' / = 0. By requiring that {kp) K \^ u is doubletraceless as any trace term would serve 
merely as a redefinition of kinematic terms and hence a field redefinition, (kp)^^ has the symmetry of Riemann 
tensor. Then we can define the decomposition of {kp)^^ in terms of its spatial and time indices, i.e., 

(k DE y k ee -2(k F ) 0j0k , (k HB f ee l(k F y kmn e^ k e lmn , 

(k DB y k ee -(k HE ) ki = \{k F )^ mn e kmn , (39) 
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with Latin indices run from 1 to 3. Moreover, we define 

as = ^tr(fc_D£;), a B = ^tr(fcjjs), (40) 

and double tracelessness gives tr(k BB + ^d_e) = 0, i.e., a = a B = —cub- So we can extract the trace term to define 

{PE) 3k = (k DE y k - aV k , -{P B ) 3k = (k HB y k + aV k . (41) 

By using Bianchi identity (/cf)k[a^] = 0, we have tr(fcDs) = 0. With these definitions, we can rewrite the Lagrangian 
© as 

photon = \(E 2 - B 2 ) + ±a(E 2 + B 2 ) + ~ ((/3 E y k EiE k + (f3 B y k B^B k + (k DB y k E^B k ) 

+ k\ F A-B-4>k A F-B + k AF -(AxE). (42) 
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